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We report explicitly a novel family of exact PT-symmetric solitons and further study their spon- 
taneous PT symmetry breaking, stabilities and collisions in Bose-Einstein condensates trapped 
in a PT-symmetric harmonic trap and a Hermite-Gaussian gain/loss potential. We observe the 
significant effects of mean-field interaction by modifying the threshold point of spontaneous PT 
symmetry breaking in Bose-Einstein condensates. Our scenario provides a promising approach to 
study PT-related universal behaviors in non-Hermitian quantum system based on the manipulation 
of gain/loss potential in Bose-Einstein condensates. 
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Parity (P) and time-reversal (T) symmetries are fun- 
damental notions in physics. There has already been 
much attention to these systems which admire the com- 
bined PT symmetry but do not obey V and T symme- 
tries separately. Despite the fact that PT-Hamiltonian 
can, in general, be non-Hermitian, the pioneering work 
of Bender and Boettcher [l[ showed that a family of 
non-Hermitian Hamiltonian with VT symmetry, H = 
p + x 2 (ix) e , can still show entirely real spectra whereas 
the generated dynamics is (pseudo) unitary correspond- 
ing to the unbroken VT symmetry. Their results have in- 
spired considerable scientific attention devoted to study 
of non-Hermitian PT-symmetric Hamiltonian in many 
aspects, such as PT-symmetric classical-mechanics the- 
ory 0], 'PT-symmetric quantum mechanics 0, PT- 
symmetric quantum field theory pseudo- Hermitian 
quantum mechanics Lie algebras M, complex crys- 
tal and PT-symmetric wave chaos [|| etc. Some re- 
view papers (|{|, for example) present an overview of the 
theoretical and experimental progress on the study of 
non-Hermitian PT-symmetric Hamiltonian. 

For the case considered here, given that the action of 
the parity P is linear and has the effect p — > —p and x —> 
—x whereas the time operator T is antilinear and has the 
effect p — > — p, x — > x, and i — > —i, it then follows that a 
necessary (but not sufficient) condition for a Hamiltonian 
to be VT symmetric is V(x) = V*{— x). In other words, 
VT symmetry implies that the real and imaginary parts 
of the PT-symmetric potential V(x) = Vr(x) + iVi(x) 
should be spatially symmetric and antisymmetric, re- 
spectively, i.e., Vr(x) = Vr(— x) and V\{x) = —Vi(—x). 
It turns out that the PT-symmetric condition is sufficient 
to guarantee that the energy spectrum is real and time 
evolution is unitary, the condition of Dirac Hermiticity is 
not necessary 

There have been several experimental and theoretical 
studies on one-dimensional Schrodinger equation with 



different PT-symmetric potentials such as the periodic 
potentials non-periodic potentials 11 1 and PT- 

symmetric lattices 12J. Some experiments in the frame- 
work of wave optics based on Al x Gai_z [l3| and Fe-doped 
LiNbC>3 [14[ provide experimental observation of the be- 



haviors in the both passive and active PT-symmetric op- 
tical coupled two-channel systems. More recently, Mus- 
slimani et al. fisj ] have shown that the optical nonlin- 
earity can shift the PT-symmetric threshold and in turn 
allow nonlinear eigenmodes with real eigenvalues to exist 
in the Scarff II potential [16[ and periodic potentials [17 1 . 



Very interestingly, by using the recently realized ex- 
perimental techniques [18l |. PT-symmetric gain/loss po- 
tential is very promising to be experimentally realizable 
with BECs in the near future. Moreover, a dissipative 
source (non-PT-symmetric potentials, i.e., the gain/loss 
distribution chosen Gaussian function) have been 
used to generate nonlinear coherent excitations in Bose- 
Einstein condensates (BECs) [HI- As is well known 
that, for ultra-cold alkali atoms, the confining potential is 
usually well approximated with the quadratic form [l9| . 
Thus, the natural question appeared is whether the spon- 
taneous VT symmetry breaking and stability of solitons 
are universal. The goal of this letter is to give the positive 
answer to this question by showing that PT-symmetric 
soliton-type solutions are pretty natural in BECs trapped 
by a harmonic potential and synthetic PT-symmetric 
gain/loss potential. 

In this Letter, we find a new class of exact PT- 
symmetric solitons and further study their spontaneous 
VT symmetry breaking, stabilities and collisions in 
BECs with a PT-symmetric harmonic trap and Hermite- 
Gaussian gain/loss potential. In particular, the mean- 
field interaction drastically modifies the threshold point 
of spontaneous VT symmetry breaking of solitons in 
BECs. These results provide a promising approach to- 
wards observation of other PT-related universal features 
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in new classes of PT-synthetic system. 

We concentrate on a quasi-one-dimcnsional (ID) BECs 
in the mean-field regime described by the ID Gross- 
Pitaevskii (GP) equation with complex PT-symmetric 
potential 
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here the complex amplitude distribution ^ = ^(z,t) 
is macroscopic order parameter, the nonlinearity ^id = 
2a s hu± is the effective ID coupling strength, where a s 
stands for the s-wave scattering length which can be 
changed through Feshbach resonance, and wj_ denotes 
the transverse confining frequency The external poten- 
tial is chosen as a harmonic trap Vr(z) = (l/2)muj 2 z 2 
with m being the atomic mass and lo z is the longitu- 
dinal confining frequency (l9| . The gain/loss potential 
Vi(z) is phcnomenologically incorporated to account for 
the manipulation of gain/loss potential in BECs by a 
novel experimental microscopy technique [l8j . 

After the dimensionless transformations of the den- 
sity, length, time and energy measured in units of 
g/(2a s ), aj_ = y / h/(muj±), ojJ_ , and huj±, we arrive at 
an effective ID GP equation with non-Hermitian PT- 
symmetric potential in the dimensionless variables 
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which is associated with SC/Sip* = in which the La- 
grangian density can be written as C = ^(V'V't ~ 4'*' l Pt) + 
|V>*I 2 + 2[w R (z) + wi(z)]|t/>| 2 + 5 |V| 4 , where v K (z) = 
(l/2)w 2 z 2 with the trap parameter ui = uj z /uj±, v\(z) = 
Huj±_Vi(z), and g is a interaction parameter. Based on the 
previous discussion, the real and the imaginary compo- 
nents of the 'PT-symmetric potential should satisfy the 
relations: vr(z) = vr(— z) and Vj(— z) = —vj(—z), re- 
spectively. We here concentrate on the attractive inter- 
action g < 0. The repulsive interaction g > does not 
pose new changes and will be considered elsewhere. 

We here focus on spatially localized soliton solutions 
of Eq. ^ for which \im^ z ^ 00 ip(z,t) = 0. Our goal is 
to seek a class of soliton-type stationary solutions in the 
form ?p(z,t) = 4>{z) cxp(— ifit) with fi being the chemical 
potential and the complex function 4>(z) € C[z] obeying 
the stationary GP equation with PT-symmetric poten- 
tial 
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Let us now consider the more interesting and physi- 
cally relevant PT-symmetric potential with the external 
trapping potential wr(z) and the gain/loss potential v\(z) 
being of the harmonic trap v-r(z) = (l/2)ui 2 z 2 and a fam- 
ily of the Hcrmite-Gaussian distributions 

i*JZ 2 

v\,n{z)=wuju} [2nH^i( v / w2;)-v^2i? n (v / wz)]e *~, (4) 



where wo is a constant parameter which can be con- 
trolled experimentally, the even integer index n (i.e., 
n = 0,2,4,...) is necessary to make sure v\^ n {z) be the 
odd function for the PT-symmetric potential and deter- 
mines the number of nodes across the Gaussian envelope, 
H n {z) = (-l) n e z2 d n e- z2 /dz n stands for the Hermite 
polynomial. For the odd integer index n, the function 
vi,n(z) given by Eq. (j4]) does not obey the PT-symmetric 
condition, so we do not considered these cases here. 

As a consequence, we find that there exists a family of 
exact PT-symmetric soliton-type solutions as 



V2w Q 



H n (y/ujz)e exp[itp n (z)], (5) 



where n = 0, 2, 4, g < 0, the chemical potential satis- 
fies the condition fi = uj\n + 1/2], and the phase ip n (z) is 
given by tp n (z) = (2/3)w Q f*[H n (y/ujs) cxp(~ uj s 2 /2)}ds. 
It is easy to show that these solutions is well localized 



since 4> n (z)\ 



0. 



In the following, without loss of generality, we will 
study the cases of n = and n = 2, respectively, for their 
spontaneous PT symmetry breaking, stabilities and col- 
lisions properties. Thus, we can write the gain/loss terms 
in the forms of 

vi,o(z) = — wowzexp (— ujz 2 /2) , (6a) 
vi >2 (z) = w u(10z - Aujz 3 ) exp (-wz 2 /2) , (6b) 

for n — and n ~ 2, respectively The corresponding 
exact PT-symmetric solitons should be written as 
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(2uiz 2 - l)e ^ exp[i(p 2 {z)}, (7b) 



for n ~ and n = 2, respectively, where 
erf (a;) denotes error function and ^2(2) = 
(2/3)w [v / 27r/o;erf(^V2z) - 4zexp(-wz 2 /2)]. 

As is well known that, for our studied PT-symmetric 
potentials, the linear spectrum can be entirely real as 
long as the wq is operated below the phase transition 
point (unbroken PT symmetry) for a fixed value of u>. 
Above this so-called PT-symmetric threshold, a phase 
transition occurs and the linear spectrum enters the com- 
plex domain. 

Notice that in our studied system, the correspond- 
ing linear spectrum problem associated with the PT- 
symmetric potential (|6a|) . exhibits an entirely real spec- 
trum provided that, 



Kl < 1-945, 



(8) 



for the case of 5 = — 1, u = 1 and n = 0. A single-peak 
soliton-type stationary state (|7a[) in the PT-symmetric 
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FIG. 1: (color online). Intensity evolution of a single-peak 
soliton state in a PT-symmetric potential (|6a|l for g = — 1 
and cj = 1. (a) iuo = 0.8. (b) iuo = 1.0. The inset depicts 
the real (top) and imaginary (bottom) component of "PT- 
symmetric potential. 



potential (pa|) for the case of g = —1, u = 1 and 
Wo = 0.8 is shown in Fig. [TJa). The stability of this 
PT-symmetric soliton has been confirmed by numerical 
simulation of Eq. ([2]) using beam propagation methods 
with adding random noise on both amplitude and phase 
of exact PT-symmetric soliton state. Thus for this given 
wo, one can see that PT-symmetric soliton states can 
be found with real energy in this system. Very interest- 
ingly, one finds that even if the potential (|6ap has be- 
low the phase transition point (wq = 1.0, its single par- 
ticle spectrum is real), single-peak soliton has already 
been unstable with breaking PT symmetry and imagi- 
nary energy. Fig. [ljb) shows the numerical simulation 
of Eq. © for the case of g = — 1, a; = 1 and Wq = 1.0 
by using the same method as in Fig. HJa). This means 
that, the single-peak soliton itself can modify the real 
part of potential through the mean field interaction term 
in Eq. @. Thus, this new effective trapping potential 
shifts the PT threshold wo of Eq. © towards a smaller 
value and in turn damages this stationary soliton state 
against breaking PT symmetry. For the obtained single- 
peak PT-symmetric soliton (|7ap we find that the quan- 
tity S = (i/2)(<M>l - (t> x <t>*) = 4 W 3/(9| ff |)exp(-3wz 2 /2) 
associated with the transverse power flow density is pos- 
itive with luq > 0, which means that the power always 
flows in one direction, i.e., from the gain to loss region. 

We next investigate multi-peak soliton-type stationary 
state ([7b| and their dynamics supported with the PT- 
symmetric potential (|6b[) . The linear spectrum proper- 
ties of such a PT-symmctric potential can be understood 
by examining the corresponding linear spectrum problem 
of Eq. ((3]), And one finds that purely real spectrums are 




FIG. 2: (color online). Intensity evolution of a multi-peak 
soliton state in a PT-symmetric potential (|6b[l for g = — 1 
and u) — 1. (a) wo — 0.1. (b) wo = 0.2. The inset depicts the 
imaginary component of PT- symmetric potential. 

possible in the range 

\w \ < 0.00002, (9) 

for the case of g = — 1, lj = 1 and n = 2. Since 
the similar idea holds for such PT-symmctric poten- 
tial (|6b[) , by numerical simulation of Eq. ^ using beam 
propagation methods, however, one finds a stable PT- 
symmetric multi-peak soliton-type state even for the pa- 
rameters of g = — 1, u> = 1 and wq = 0.1 as is shown 
in Fig. [2ja). Thus, very interestingly, in contrast to the 
stable PT-symmetric single-peak soliton, here, the multi- 
peak soliton itself eventually, shifts the PT threshold 
wo of Eq. ([9]) towards a rather larger value and makes 
soliton-type state preserve PT symmetry very well by 
drastically modifying the real part of potential through 
the mean field interaction term. With increasing the 
wq further, the multi-peak soliton will be unstable with 
breaking PT symmetry and imaginary energy. Fig. [5{b) 
shows the numerical simulation of Eq. ([2]) for the case of 
g = — 1, u> = 1 and wo = 0.2 by using the same method 
as former cases. 

Finally, we need to check if this PT-symmetric soliton 
indeed behave as solitons. To this end, we take initial 
conditions in the form of a superposition of two single- 
peak (|7aj) or multi-peak (|7b|) solitons with symmetric 
central positions. In such a case, and given that the 
single-peak or multi-peak solitons were found above to 
be robust objects behaving similarly to solitons of an in- 
tegrable system, one may expect that the solitons would 
perform harmonic oscillations in the presence of the suf- 
ficiently large parabolic trap even with PT-symmetric 
potentials. As shown in the Fig. Eta) for the case of 
g = —1, cu = 1 and wq = 0.3, two single-peak solitons 
propagate in opposite directions and undergo a pretty 
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in BECs with TT-symmetric gain/loss potential. Their 
spontaneous VT symmetry breaking, stability, as well as 
collision dynamics of such T-'T-symmetric solitons were 
examined in detail. We observe the significant effects of 
mean-field potential by modifying the threshold point of 
spontaneous VT symmetry breaking. Our results may 
intrigue a new class of 'PT-synthctic universal proper- 
ties study that relies on manipulations of macroscopic 
matter-wave field with BECs. 

This work was supported by the NSFC under grants 
Nos. 60821001/F02, 10874235, 10934010, 60978019, 
the NKBRSFC under grants Nos. 2009CB930701, 
2010CB922904 and 2011CB921500. 



FIG. 3: (color online). Dynamic evolution of collisions of 
two single-peak and multi-peak solitons in a PT-symmetric 
potentials (|6a[l and (|6b|) for g = — 1 and w = 1. Here (a) 
Wo = 0.3 and (b) wq — 0.8 for two single-peak solitons. (c) 
Wq = 0.2 and (d) u>o = 0.25 for two multi-peak solitons. 



clastic collision and the two solitons remain unscathed 
after the collision. The stability of this superposition of 
two single-peak solitons has been confirmed by numerical 
simulation of Eq. ([2]) using beam propagation methods as 
before. But when the gain/loss parameter wq becomes 
larger (for example, w = 0.8), two single-peak solitons 
undergo a inelastic collision as shown in the Fig.^b). A 
superposition of two multi-peak solitons (|7b|) with sym- 
metric central positions is also found to be robust objects 
behaving similarly to above and undergos a elastic col- 
lision as shown in the Fig. [3^c) for the case of g = —1, 
u = 1 and wq = 0.2. But when the gain/loss parameter 
wq becomes a little larger (for example, wq = 0.25), two 
multi-peak solitons become unstable clearly as shown in 
the Fig. [3jd). All these findings certainly indicate that 
many kinds of robust stable nonlinear modes can be sup- 
ported by PT-symmetric potential in BECs. 

Experimental realization of the ideas presented here 
is promising. A condensate in a dilute vapor of sodium 
or rubidium atoms is placed in a cigar-shaped potential 
may provide an accessible laboratory to experimentally 
observe spontaneous VT symmetry breaking and stabil- 
ity of solitons. Such an experiment would require precise 
control over the rates of gain/loss process in the con- 
densate where atoms are injected continuously into the 
condensate in left side, whereas they are removed simul- 
taneously from the condensate in right side. The com- 
bination of high spatial resolution with in situ detection 
[l8l | will hopefully provide new possibilities for the prepa- 
ration, the manipulation and the characterization of non- 
Hcrmitian VT symmetry condensate. 

In summary, we have demonstrated that a novel class 
of exact PT-symmctric solitons can be well supported 
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